home *** CD-ROM | disk | FTP | other *** search
/ Shareware Overload Trio 2 / Shareware Overload Trio Volume 2 (Chestnut CD-ROM).ISO / dir24 / aprs501a.zip / MAPFIX.BAS < prev    next >
BASIC Source File  |  1994-05-12  |  55KB  |  1,194 lines
  1. REM MAPFIX.bas PROGRAM.  SEE EXPLAINATION BELOW
  2. REM
  3. Ver$ = "4.05"
  4. MaxNumMAPS = 120'was 99 Current maximum number of maps loaded by APRS
  5. MaxNumPoints = 1500 'was 1000
  6. MaxNumLABELS = 99 'was 99
  7. MaxNumLines = 900
  8. REM $DYNAMIC
  9. GOTO BEGIN
  10.  
  11. Info: COLOR 15, 4: CLS
  12.    PRINT " MAPFIX.bas VERSION "; Ver$; " PROGRAM FOR FIXING APRS MAPS": PRINT
  13.    PRINT " This program has evolved significantly.  305 added a TRIM command to eliminate"
  14.    PRINT " all points outside an area.  3.07B adds W7KKE's GPS track history overlays and"
  15.    PRINT " digitizer code for using a drawing tablet.  308 fixed a filename bug.  Now in"
  16.    PRINT " 312, it will take USGS data from the 2,000,000:1 CD ROM and build APRS maps!"
  17.    PRINT " BUT NOT EASILY!  For a 50 mile map, there are 5 MBytes of USGS data which must"
  18.    PRINT " be filtered to get down to the 10K APRS size for the same area! The final steps"
  19.    PRINT " are all manual and take almost as long as developing an APRS map from scratch!"
  20.    PRINT
  21.    PRINT " CAUTION, THIS PROGRAM IS NOT PERFECT... KEEP BACKUPS!  Do a little at a time!"
  22.    PRINT
  23.    PRINT " Although MAPFIX has commands to make map modifications easier, a text EDITOR is"
  24.    PRINT " still useful for whole scale rearranging of points and features in a map file. "
  25.    PRINT
  26.    PRINT " MAPFIX uses two cursors.  The normal yellow APRS cursor, and a White MapPoint"
  27.    PRINT " which will be the next point to be processed.  ALT Keys allow you to MOVE the"
  28.    PRINT " MapPoint to the cursor, ADD a new point at the cursor, or DELETE the MapPoint."
  29.    PRINT " MAPFIX.bas shows you Deg/Min, Decimal, and APRS values of the cursor position."
  30.    PRINT
  31.    PRINT " ALSO NOTE THAT THE LIMITS IN APRS ARE 1000 POINTS, 99 FEATURES, and 99 LABELS!"
  32.    PRINT " If you need more points, features or labels, begin another map.";
  33.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  34.    GOSUB GetChar: SOUND 800, 3
  35.   
  36. Info2:   COLOR 15, 8: CLS
  37.    PRINT " PAGE 2 INSTRUCTIONS:  More about new features in version "; Ver$
  38.    PRINT
  39.    PRINT " With this new MAPFIX.bas, you can not only modify features by moving, adding"
  40.    PRINT " or deleting POINTS, but you can now add and kill FEATURES too, ie: roads,"
  41.    PRINT " rivers, borders, etc, from within the program.  In addition to the new KILL"
  42.    PRINT " and TRIM commands, the G key will move the cursor (GOTO) to the MapPoint,"
  43.    PRINT " and the F key will FIND the MapPoint nearest the cursor location.  If the"
  44.    PRINT " MapPointer and FeatureName get out of sequence, the RESET command may fix"
  45.    PRINT " them, but you should save the file immediately and check it with an editor."
  46.    PRINT
  47.    PRINT " I find the capability to delete points very useful when making larger area"
  48.    PRINT " maps from several smaller detail maps.  First, I run MAPCNVRT.bas to convert"
  49.    PRINT " all of the smaller maps to new temporary files with the new origin of the new"
  50.    PRINT " larger map.  Then I use the KILL command in MAPFIX to eliminate all minor "
  51.    PRINT " roads, features and labels and then the DELETE POINT command to remove all"
  52.    PRINT " inconsequential minor points from the roads that will not be needed"
  53.    PRINT " at the larger scale.  Then I use the editor to combine all of the points and"
  54.    PRINT " labels into the new file."
  55.    PRINT
  56.    PRINT " A new MAPLIST command shows your MAPLIST.map file; and the OTHER MAPS command"
  57.    PRINT " shows all MAP borders so you can see how your new map fits in.  You may use F3"
  58.    PRINT " and F4 keys to select smaller or larger map borders to draw."
  59.    PRINT
  60.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  61.    GOSUB GetChar: SOUND 800, 3
  62.   
  63. Info3:   COLOR 15, 3: CLS
  64.    PRINT " PAGE 3 INSTRUCTIONS:  Using GPS Track History Files to draw maps!"
  65.    PRINT
  66.    PRINT " To aid in creating accurate maps, W7KKE in California added routines to MAPFIX"
  67.    PRINT " so that you can overlay a Track History file onto the map you are constructing."
  68.    PRINT " This is an excellent tool for correcting your maps to real GPS data."
  69.    PRINT ""
  70.    PRINT " After you have loaded your map, type alt-G (GPS Track History') and enter the"
  71.    PRINT " history filename.  This will overlay the track history file.  You may then"
  72.    PRINT " use the normal MAPFIX.bas routines to move map segments and add so that the"
  73.    PRINT " map will agree with the GPS data contained in the track history file.  This"
  74.    PRINT " is especially useful with the history files saved by a laptop during mobile"
  75.    PRINT " GPS operations."
  76.    PRINT
  77.    PRINT " CAUTION:  Since GPS data is only accurate to 100 yards due to the effects of"
  78.    PRINT " Selective Availability, I would avoid using GPS data explicitely below about"
  79.    PRINT " the 2 mile range.  For this reason, I make the size of the GPS positions "
  80.    PRINT " expand below the 2 mile range to roughly approximate the size of the 100 yard"
  81.    PRINT " error circle."
  82.    PRINT
  83.    PRINT " Also note that you can now START a NEW map from scratch, without using the"
  84.    PRINT " text editor to enter the initial configuration data.  Just type NEW instread "
  85.    PRINT " of a MAPfilename when starting up the program."
  86.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  87.    GOSUB GetChar: SOUND 800, 3
  88.  
  89.  
  90. Info4:   COLOR 15, 7: CLS
  91.    PRINT " PAGE 4 INSTRUCTIONS FOR USING A DIGITIZER:": PRINT
  92.    PRINT " In version 3.07B, MAPFIX.bas can now interface to a digitizer tablet or table"
  93.    PRINT " so that maps can be drawn directly from paper maps.  These routines were"
  94.    PRINT " developed by W7KKE and added to MAPFIX.bas in December 1993."
  95.    PRINT
  96.    PRINT " To use a digitizer, first you must hit the ALT-O command to open the COMM"
  97.    PRINT " PORT for the digitizer.  This command also lets you test the digitizer while"
  98.    PRINT " testing the alignment of the map on the digitizer surface.  It then prompts"
  99.    PRINT " you to identify the upper left and lower right corners of the map, in order"
  100.    PRINT " to calibrate the digitizer to the latitude, longitude and scale of the map."
  101.    PRINT
  102.    PRINT " From this point on, the button on the digitizer mouse is almost identical to"
  103.    PRINT " the ALT-A command for ADDing a point.  To start a new map feature, however,"
  104.    PRINT " for the digitizer, you DO NOT use the ALT-N NEW command, but you should use"
  105.    PRINT " the ALT-B BEGIN command.  For More information, see the README.DIG file."
  106.    PRINT
  107.    PRINT " To speed up the map drawing during editing, I no longer erase and re-draw"
  108.    PRINT " the entire map with each new point.  I simply draw just the new line segment."
  109.    PRINT " Sometimes, especially when you move, or add a line, this leaves an old line"
  110.    PRINT " segment, where there actually is no longer a line.  You can always celan up"
  111.    PRINT " the map by just hitting the space bar to force a new map..."
  112.    PRINT
  113.  
  114.   
  115.    Display$ = "UNKnown"
  116.    RETURN
  117.  
  118. GetChar: a$ = "": DO UNTIL a$ <> "": a$ = INKEY$: LOOP: RETURN
  119.  
  120. BEGIN: GOSUB Info:
  121.    PRINT " HIT ANY KEY to proceed onto the HELP screen...";
  122.    GOSUB GetChar
  123.  
  124.    DIM x%(5 * MaxNumPoints), y%(5 * MaxNumPoints)
  125.    REM MAP coordinates               **** THESE ARE BIGGER THAN APRS ***
  126.    DIM LN$(MaxNumLines) ' (no limit in APRS)  **** SO YOU CAN MANIPULATE BIG MAPS
  127.    nn = 2 * MaxNumLABELS
  128.    DIM ML$(nn), MLa(nn), MLo(nn), MLr(nn) 'Map Labels, lengths and coordinates
  129.    nn = 2 * MaxNumMAPS
  130.    DIM MapName$(nn), LATcen(nn), LONcen(nn), MapMax(nn), Comment$(nn)
  131.    RdsOn = -1: Labls = -1: Tags = -1: KP = 1: Changed = 0: MapSize = 256
  132.    i = 1000
  133.    DIM HLAT(i), HLONG(i)'For lat/longs from big GPS history files
  134.  
  135. INIT: ON ERROR GOTO ErrorTrap
  136.    ScrnType$ = "EGA": Ycen = 200: Yfactr = 1: YfacTXT = 350 / 350: SCREEN 9
  137.    IF ScrnType$ = "EGA" THEN COLOR 15, 0
  138.    REM ScrnType$ = "CGA": Yfactr=200/400:Ycen = 200*Yfactr: SCREEN 2
  139.    ReDraw = -1
  140.   
  141. Display$ = "HELP": GOSUB HELP: GOSUB LoadMap
  142. REM ON ERROR GOTO 0
  143.  
  144. Main: GOSUB DrwMPaCur
  145.   DO
  146. GoAgain: Fault = 0
  147.     IF Digitizer THEN
  148.          IF LOC(1) > 9 THEN
  149.             GOSUB GetXY: GOSUB Cursor
  150.             IF Btn <> 3 THEN GOSUB AddPoint
  151.          END IF
  152.     END IF
  153.     a$ = INKEY$
  154.     IF a$ <> "" THEN
  155.         a$ = UCASE$(a$): Key$ = a$
  156.         IF a$ = "S" THEN GOSUB labels
  157.         IF a$ = "L" THEN Labls = NOT Labls
  158.         IF a$ = "T" THEN Tags = NOT Tags
  159.         IF a$ = "F" THEN GOSUB FindPoint
  160.         IF a$ = "G" THEN GOSUB CurToPoint: GOSUB CurDrwMap
  161.         IF a$ = "H" THEN
  162.            IF Display$ <> "HELP" THEN
  163.               GOSUB HELP
  164.            ELSE GOSUB Info
  165.               LOCATE 25, 1: PRINT " H for HELP or SPACE BAR for map..."; : a$ = ""
  166.            END IF
  167.         END IF
  168.         IF a$ = "B" THEN GOSUB BoxPPD
  169.         IF a$ = "U" THEN GOSUB GetUSGS
  170.         IF a$ = "D" THEN GOSUB MapDIR
  171.         IF a$ = "M" THEN GOSUB ListMAPlist
  172.         IF a$ = "O" THEN GOSUB DrwAndShow
  173.         IF a$ = "N" THEN GOSUB NextLine: GOSUB Cursor
  174.         IF a$ = "P" THEN GOSUB Previous: GOSUB Cursor
  175.         IF a$ = "Q" THEN GOSUB QUIT
  176.         IF a$ = "R" THEN Z = 2: LNptr = 1
  177.         IF a$ = "T" THEN GOSUB Scrunch
  178.         IF a$ = " " THEN Display$ = "MAP": ReDraw = -1: USGS = 0: GOSUB DrwMPaCur
  179.         IF a$ = "+" THEN Z = Z + 1: GOSUB MapPoint ' moves to next map point
  180.         IF a$ = "-" THEN Z = Z - 1: GOSUB MapPoint ' moves backwards
  181.         IF a$ = CHR$(18) THEN ReDraw = NOT ReDraw: GOSUB ReDraw
  182.        
  183.         B$ = "": IF LEN(a$) = 2 THEN B$ = RIGHT$(a$, 1): REM process arrow & special keys
  184.         IF B$ = "I" THEN RS = RS * 2: GOSUB CurDrwMap: REM change scale
  185.         IF B$ = "Q" THEN RS = RS / 2: GOSUB CurDrwMap
  186.         IF B$ = CHR$(132) THEN RS = RS * 8: GOSUB CurDrwMap: REM change scale by factor of 4
  187.         IF B$ = "V" THEN RS = RS / 8: GOSUB CurDrwMap
  188.         IF B$ = "G" THEN GOSUB CurDrwMap 'Home key
  189.         IF a$ = "7" THEN CDX = LONo: CDY = LATo: GOSUB DrwMPaCur 'ShiftHOME
  190.         IF B$ = "O" THEN CDX = LONcen: CDY = LATcen: GOSUB DrwMPaCur 'End Key
  191.         IF B$ = "M" THEN CPX = CPX - 4 / (Sfac): GOSUB Cursor
  192.         IF B$ = "K" THEN CPX = CPX + 4 / (Sfac): GOSUB Cursor
  193.         IF B$ = "H" THEN CPY = CPY + 4 / (Sfac): GOSUB Cursor
  194.         IF B$ = "P" THEN CPY = CPY - 4 / (Sfac): GOSUB Cursor
  195.         REM Here are the special MapFIx routines
  196.         IF B$ = CHR$(30) THEN GOSUB AddPoint               'alt-ADD point
  197.         IF B$ = CHR$(48) AND Digitizer THEN GOSUB NewFeature'alt-BEGIN
  198.         IF B$ = CHR$(34) THEN GOSUB LoadHst                'alt-GPS hstry file
  199.         IF B$ = CHR$(50) THEN GOSUB MakePT: IF ReDraw THEN GOSUB DrawMap 'MOVE point to cursor
  200.         IF B$ = CHR$(32) THEN GOSUB DelPT                  'alt-DELete point
  201.         IF B$ = CHR$(38) THEN GOSUB AddLabel               'alt-ADD LABEL
  202.         IF B$ = CHR$(46) THEN GOSUB NewCenter              'alt-CENTER
  203.         IF B$ = CHR$(36) THEN GOSUB Join                   'alt-JOIN
  204.         IF B$ = CHR$(37) THEN GOSUB KillF                  'alt-KILL Feature
  205.         IF B$ = CHR$(19) THEN GOSUB MapRange               'alt-RANGE
  206.         IF B$ = CHR$(20) THEN GOSUB TRIM                   'alt-TRIM
  207.         IF B$ = CHR$(49) THEN GOSUB NewFeature             'alt-NEW Feature
  208.         IF B$ = CHR$(24) THEN GOSUB DigiInit: GOSUB DrawMap'alt-OPEN dgtzr COM
  209.         IF B$ = CHR$(31) THEN GOSUB Scrunch                'alt-SCRUNCH
  210.         IF B$ = CHR$(22) THEN GOSUB GetUSGS                'alt-U
  211.         IF B$ = CHR$(61) THEN                              'F3 for smaller Maps
  212.            MapSize = MapSize / 2: IF MapSize < 1 THEN MapSize = 1
  213.            GOSUB ShowMaps
  214.            END IF
  215.         IF B$ = CHR$(62) THEN                              'F4 for larger Maps
  216.            MapSize = MapSize * 2: IF MapSize > 1000 THEN MapSize = 1000
  217.            GOSUB DrwAndShow
  218.         END IF
  219.         IF a$ = CHR$(19) THEN GOSUB SaveMap
  220.         IF a$ = CHR$(3) THEN GOSUB ChgColr
  221.            
  222.         IF a$ = "6" THEN CPX = CPX - 20 / (Sfac): GOSUB Cursor'SHIFT Cursor by 4
  223.         IF a$ = "4" THEN CPX = CPX + 20 / (Sfac): GOSUB Cursor
  224.         IF a$ = "8" THEN CPY = CPY + 20 / (Sfac): GOSUB Cursor
  225.         IF a$ = "2" THEN CPY = CPY - 20 / (Sfac): GOSUB Cursor
  226.         
  227.     END IF
  228.   LOOP
  229.   SYSTEM 'you should never get here
  230.  
  231. ReDraw: LOCATE 1, 30
  232.         IF ReDraw THEN PRINT "REDRAW ENABLED":  ELSE PRINT "NO ReDraw...  "
  233.         RETURN
  234.  
  235. QUIT: a$ = "Y"
  236.       IF Changed THEN
  237.          GOSUB BoxLine23
  238.          PRINT "**** MAP HAS BEEN MODIFIED"; Changed; "TIMES BUT NOT SAVED!!!  SAVE NOW? (Y)";
  239.          INPUT a$
  240.       IF UCASE$(a$) <> "N" THEN GOSUB SaveMap
  241.       END IF
  242.       SYSTEM
  243.  
  244. TRIM: GOSUB BoxLine23
  245.       CLS : PRINT "TRIM ALL POINTS AND LABELS OUTSIDE OF MAPRANGE"
  246.       PRINT
  247.       PRINT "This command will remove all points and labels that are outside of the white"
  248.       PRINT "map border.  You can change the location of this map border by using"
  249.       PRINT "the CENTER command (alt-C) and by changing the RANGE using alt-R."
  250.       PRINT : PRINT
  251.       PRINT "No map feature will be completely eliminated..."
  252.       PRINT
  253.       PRINT "The first and last point of any FEATURE will be retained, so the"
  254.       PRINT "result will be long single lines for all FEATURES outside the map border."
  255.       PRINT "Use the KILL FEATURE (alt-K) to eliminate those lines and use the MOVE"
  256.       PRINT "command (alt-M) to move any far away points closer to the border."
  257.       PRINT : PRINT
  258.       PRINT "You might consider stopping now and doing a SAVE (ctrl-S) before proceeding."
  259.       PRINT
  260.       PRINT "ALSO, THIS DOES NOT WORK FOR POINTS WITH NEGATIVE VALUES!  Be sure  your"
  261.       PRINT "selected area is below and to right of ORIGIN.  If not, run MAPCNVRT.bas."
  262.       PRINT : PRINT
  263.       INPUT "Are you ready to proceed? (Y/N) (N)"; ans$
  264.       GOSUB DrawMap
  265.       IF UCASE$(ans$) <> "Y" THEN RETURN
  266.       C = 0: LOCATE 23, 1: PRINT "Processing...";
  267.       REM dx and dy are num pix of center of map
  268.       REM bx and by are borders of map based on MapRng
  269.       by = ppdV * MapRng / 60
  270.       bx = by / Lfac
  271.       FOR Z = 1 TO nmp - 4
  272.          IF x%(Z) = 0 THEN Z = Z + 2
  273.          IF x%(Z) > dx + bx OR y%(Z) > dy + by THEN bad = 1 ELSE bad = 0
  274.          IF x%(Z) < dx - bx OR y%(Z) < dy - by THEN bad = 1
  275.          IF bad AND x%(Z - 1) <> 0 AND x%(Z + 1) <> 0 THEN
  276.             GOSUB DelPT: Z = Z - 1
  277.             C = C + 1
  278.          END IF
  279.       NEXT Z
  280.       LOCATE 23, 1: PRINT "Now removing labels...";
  281.       FOR i = 1 TO nml: REM now eliminate all labels outside
  282.           bad = 0: Xm = MapRng / (60 * Lfac): Ym = MapRng / 60
  283.           IF MLo(i) > LONcen + Xm OR MLa(i) > LATcen + Ym THEN bad = 1
  284.           IF MLo(i) < LONcen - Xm OR MLa(i) < LATcen - Ym THEN bad = 1
  285.           IF bad = 1 THEN
  286.              FOR j = i TO nml
  287.                  ML$(j) = ML$(j + 1): MLa(j) = MLa(j + 1)
  288.                  MLo(j) = MLo(j + 1): MLr(j) = MLr(j + 1)
  289.              NEXT j: nml = nml - 1: PRINT ".";
  290.           END IF
  291.       NEXT i
  292.       GOTO DrawMap
  293.  
  294.  
  295. FindPoint: CurX = INT(.5 + dx + (CUX - 320) / KP)
  296.            CurY = INT(.5 + dy + (CUY - Ycen) / KP)
  297.      GOSUB BoxLine23: PRINT "SEARCHING THROUGH ALL POINTS IN FILE...";
  298.      Pass = 1: SaveZ = Z: SaveLNptr = LNptr
  299. Agn: FOR j = 0 TO 30            ' Go through abt 20 times lookin pt.
  300.         IF j > 10 THEN j = j + 1' first with 0 delta, then bigger
  301.         LNctr = LNptr
  302.         PRINT ".";
  303.         FOR i = Z + 1 TO nmp
  304.             IF x%(i) = 0 THEN LNctr = LNctr + 1
  305.             IF x%(i) > CurX - j AND x%(i) < CurX + j THEN
  306.                IF y%(i) > CurY - j AND y%(i) < CurY + j THEN
  307.                   Z = i: LNptr = LNctr: GOSUB CurToPoint
  308.                   j = 99: i = nmp
  309.                END IF
  310.             END IF
  311.         NEXT i:
  312.     NEXT j
  313.     IF j < 99 AND Pass = 1 THEN Z = 2: LNptr = 1: Pass = 0: GOTO Agn
  314.     IF j < 99 THEN PRINT "None found!": LNptr = SaveLNptr: Z = SaveZ: RETURN
  315.     GOSUB MapPoint: SavClr = 0: RETURN
  316.                   
  317. NewFeature: LOCATE 24, 1: PRINT SPACE$(27); : GOSUB BoxLine23
  318.             INPUT "Enter reference name for new feature"; a$
  319.             IF a$ = "" THEN RETURN
  320.             GOSUB Rainbow: IF abort THEN RETURN
  321.             GOSUB BeginF
  322.             GOSUB BoxLine23: LOCATE 25, 1: PRINT SPACE$(80); : LOCATE 25, 1
  323.             IF RIGHT$(Key$, 1) = CHR$(48) THEN
  324.                PRINT "NOW USE DIGITIZER TO ADD NEW POINTS TO THIS FEATURE...";
  325.                GOSUB GetXY: GOSUB Cursor
  326.             ELSE
  327.                PRINT "NOW MOVE CURSOR AND USE ALT-A TO ADD POINTS TO THIS NEW FEATURE...";
  328.             END IF
  329.             GOSUB MakePT
  330.             RETURN
  331.  
  332. Rainbow: LOCATE 25, 1
  333.          FOR i = 0 TO 14
  334.              PRINT RIGHT$(" " + MID$(STR$(i + 1), 2), 2); "   ";
  335.              LINE (16 + i * 40, 335 * YfacTXT)-(40 + i * 40, 349 * YfacTXT), i + 1, BF
  336.          NEXT i
  337.          GOSUB BoxLine23
  338.          INPUT "Select color (4,7,10-Hwys 11-Water 12-Hwy 13-Spcl 14-City)"; B$
  339.          SavClr = VAL(B$): IF SavClr > 15 OR SavClr < 1 THEN abort = -1 ELSE abort = 0
  340.          RETURN
  341.  
  342.            
  343. BeginF: x%(nmp) = 0: y%(nmp) = SavClr   'Store feature color 0,c
  344.         LN$(LNi + 1) = LN$(LNi): LNptr = LNi'Bump up present LN$ comment
  345.         LN$(LNi) = a$: LNi = LNi + 1'Store feature name
  346.         nmp = nmp + 1: Z = nmp
  347.         nmp = nmp + 1: x%(nmp) = 0: y%(nmp) = 0'nmp points to ending 0,0
  348.         RETURN
  349.  
  350. CanclF: nmp = nmp - 2: Z = Kz
  351.         LNi = LNi - 1: LN$(LNi) = LN$(LNi + 1): RETURN
  352.  
  353. NewCenter: LATcen = CPY: LONcen = CPX: Changed = Changed + 1: GOTO CurDrwMap
  354.  
  355. MapRange: GOSUB BoxLine23: INPUT "Enter map range"; a$
  356.           IF VAL(a$) <> 0 THEN MapRng = VAL(a$)
  357.           Changed = Changed + 1: GOTO DrwMPaCur
  358.  
  359. AddPoint: nmp = nmp + 1: Z = Z + 1
  360.           FOR i = nmp TO Z STEP -1
  361.               x%(i) = x%(i - 1): y%(i) = y%(i - 1)
  362.           NEXT
  363.           GOSUB MakePT
  364.           IF SavClr = 0 AND ReDraw THEN GOTO DrawMap
  365.           s = Z - 1: LineColor = SavClr: GOTO DP
  366.  
  367. MakePT: x%(Z) = dx + (CUX - 320) / (KP * Hfac)
  368.         y%(Z) = dy + (CUY - Ycen) / KP
  369.         Changed = Changed + 1
  370.         GOTO MapPoint
  371.  
  372. CurToPoint:
  373.      CPX = CDX - (x%(Z) - dx) / ppdV
  374.      CPY = CDY - (y%(Z) - dy) / (ppdV * Yfactr)
  375.      GOTO Cursor
  376.  
  377. DelPT: GOSUB DelZ
  378.        REM if 1st pt, it stays as 1st pt
  379.  
  380.        IF x%(Z) = 0 THEN Z = Z - 1: REM if end pt, it stays as end
  381.        IF x%(Z + 1) = 0 AND x%(Z - 1) = 0 THEN 'It is LAST point
  382.           GOSUB Kline: LNptr = LNptr - 1       'So Kill Line
  383.           GOSUB DelZ                           'And Kiil it
  384.           Z = Z - 1: GOSUB DelZ: Z = Z - 1     'Kill 0,color
  385.        END IF                                  'and -1 to end point
  386.        IF B$ = CHR$(32) AND ReDraw THEN GOSUB DrawMap ELSE GOSUB MapPoint
  387.        RETURN
  388.  
  389. DelZ: nmp = nmp - 1
  390.       FOR i = Z TO nmp
  391.           x%(i) = x%(i + 1): y%(i) = y%(i + 1)
  392.       NEXT: Changed = Changed + 1: RETURN
  393.  
  394. NextLine: IF Z >= nmp - 1 THEN Z = nmp - 1: BEEP: RETURN
  395.           DO UNTIL x%(Z) = 0: Z = Z + 1: LOOP
  396.           IF Z < nmp - 1 THEN Z = Z + 1: LNptr = LNptr + 1
  397.           SavClr = 0: GOTO MapPoint
  398. Previous: DO UNTIL Z = 1 OR x%(Z) = 0: Z = Z - 1: LOOP
  399.           IF Z > 3 THEN Z = Z - 1: LNptr = LNptr - 1
  400.           SavClr = 0: GOTO MapPoint
  401.  
  402. KillF: GOSUB Find1st: REM Stop at Beginning (0) point of the feature to kill
  403.        ni = Bi + 1' Now scan for next feature
  404.        DO UNTIL x%(ni) = 0: ni = ni + 1: LOOP
  405.        REM now move down rest of array to fill
  406.        DO UNTIL ni = nmp + 1
  407.            x%(Bi) = x%(ni): y%(Bi) = y%(ni)
  408.            Bi = Bi + 1: ni = ni + 1
  409.        LOOP
  410.        nmp = nmp - (ni - Bi): y%(nmp) = 0
  411.        GOSUB Kline
  412.        GOTO DrawMap
  413.     
  414. Find1st: Bi = Z: Changed = Changed + 1
  415.          DO UNTIL x%(Bi) = 0: Bi = Bi - 1: LOOP: Z = Bi + 1
  416.          RETURN
  417.  
  418. ChgColr: GOSUB Find1st: GOSUB Rainbow: IF abort THEN RETURN
  419.          y%(Bi) = SavClr: RETURN
  420.  
  421.  
  422. Kline: FOR i = LNptr TO LNi
  423.            LN$(i) = LN$(i + 1)
  424.            NEXT i
  425.            LNi = LNi - 1
  426.            RETURN
  427.  
  428. MapPoint:
  429.      IF Z < 2 THEN Z = 2: LNptr = 1: BEEP: SavClr = 0
  430.      IF Z > nmp - 1 THEN Z = Z - 1: BEEP: SavClr = 0
  431.      IF x%(Z) = 0 THEN
  432.         IF a$ = "-" THEN
  433.              LNptr = LNptr - 1: Z = Z - 1
  434.         ELSE LNptr = LNptr + 1: Z = Z + 1
  435.         END IF: SavClr = 0
  436.      END IF
  437.      IF LNptr < 0 THEN LNptr = 0
  438.      IF Display$ = "MAP" THEN
  439.           LOCATE 22, 1
  440.           PRINT "Fture#"; LNptr; TAB(12); LEFT$(LN$(LNptr) + "            ", 12);
  441.      END IF
  442. DrwMpPt: IF Display$ <> "MAP" THEN RETURN
  443.      CIRCLE (Xtest, Ytest), 10, 0 'Erase old circle
  444.      Xtest = 320 + KP * (x%(Z) - dx) * Hfac
  445.      Ytest = Ycen + KP * (y%(Z) - dy) * Yfactr
  446.      CIRCLE (Xtest, Ytest), 10, 15
  447.      
  448.      LOCATE 23, 1: PRINT "MapPt#"; Z;
  449.      IF Z > 999 THEN PRINT TAB(13); "val:";  ELSE PRINT TAB(12); "vals:";
  450.      PRINT TAB(17); x%(Z); TAB(23); y%(Z)
  451.      RETURN
  452.  
  453. AddLabel: nml = nml + 1
  454.           MLa(nml) = CPY: MLo(nml) = CPX
  455.           GOSUB BoxLine23: INPUT "Enter Label Name"; a$: ML$(nml) = a$
  456.           GOSUB BoxLine23: INPUT "Begin displaying label at what range?"; a$
  457.           a = VAL(a$): IF a <> 0 THEN MLr(nml) = a:  ELSE MLr(nml) = 2048
  458.           Changed = Changed + 1: GOTO labels
  459.  
  460. BoxLine23: LOCATE 23, 1: PRINT SPACE$(80); : LOCATE 23, 1: RETURN
  461.  
  462. ErrorTrap: Fault = ERR: 'Error handling routine
  463.            IF ERR = 57 THEN PRINT "  I/O-error-User-logoff"; : RESUME
  464.            IF ERR = 69 THEN PRINT "  Comm-buffer-overflow"; : RESUME
  465.            IF ERR = 53 THEN PRINT "  file-"; F$; "-not-found": CLOSE : RESUME NEXT
  466.            IF ERR = 62 THEN RESUME NEXT
  467.            IF ERR = 52 THEN RESUME NEXT
  468.            IF ERR = 55 THEN RESUME NEXT
  469.            IF ERR = 2 THEN PRINT "SYNTAX-error"
  470.            IF ERR = 70 THEN PRINT " WRITE PROTECTED!...": RESUME NEXT
  471.            IF ERR = 76 THEN PRINT "Wrong Path!": RESUME NEXT
  472.            RESET
  473.            PRINT : PRINT "Error beyond repair. Number = "; ERR;
  474.            INPUT "Hit RETURN to return to DOS"; a$
  475.            SYSTEM
  476.  
  477. MapDIR: CLS : PRINT "MAP FILES DIRECTORY": PRINT
  478.          PRINT "To display MAP files, please enter the path to your xxxxxxx.MAP files."
  479.          PRINT "For example, the default '\APRS\MAPS\*.MAP' will show all maps in the APRS"
  480.          PRINT "directory.  Similarly '*.map' will search your present QB directory."
  481.          PRINT "For any other path, enter the full file specification.": PRINT
  482.          F$ = "\aprs\MAPS\*.map"
  483.          PRINT "Enter Filespec for searching the DIRECTORY ("; F$; ")";
  484.          INPUT a$: IF a$ <> "" THEN F$ = a$
  485.          PRINT : PRINT : FILES F$
  486.          RETURN
  487.  
  488.  
  489. LoadMap: 'Maps are drawn to the default EGA resolution of 640 x 400 (350)
  490. Again: GOSUB BoxLine23
  491.        INPUT " Enter map FILENAME, or NEW, or ? for a list, or Q to quit)"; a$
  492.        a$ = UCASE$(a$): IF a$ = "" THEN GOTO Again
  493.        IF a$ = "Q" THEN SYSTEM
  494.        IF a$ = "?" THEN GOSUB MapDIR: GOTO Again
  495.        IF a$ = "NEW" THEN Key$ = "NEW": GOSUB NewMap: RETURN
  496.        a = INSTR(3, a$, "."): IF a = 0 THEN a$ = a$ + ".MAP"
  497.        MapFile$ = a$: F$ = MapFile$: OPEN F$ FOR INPUT AS #3
  498.        IF Fault = 53 THEN Fault = 0: PRINT : CLOSE #3: GOTO Again
  499.        GOSUB BoxLine23: PRINT " Loading "; F$; "..."
  500.        INPUT #3, LATo: LINE INPUT #3, LATtext$
  501.        INPUT #3, LONo: LINE INPUT #3, LONtext$
  502.        INPUT #3, ppdV: LINE INPUT #3, VS$'Pixels per degree horiz
  503.        INPUT #3, LATcen: LINE INPUT #3, LATcen$
  504.        INPUT #3, LONcen: LINE INPUT #3, LONcen$
  505.        INPUT #3, MapRng: LINE INPUT #3, MapRng$
  506.        INPUT #3, MinRng: LINE INPUT #3, MR$
  507.        LINE INPUT #3, TextLine$ ' Line of comments or instrutcitons
  508.        IF LEFT$(TextLine$, 14) = "Map generated " THEN ReDraw = 0
  509.        RS = 2 ^ INT(LOG(MapRng) / LOG(2))'Rng is intgr of VERTrng
  510.        i = 0: LNi = 0:
  511.      
  512.      DO WHILE NOT EOF(3)
  513.         i = i + 1: INPUT #3, x%(i), y: y%(i) = y * Yfactr
  514.         IF x%(i) = 0 AND NOT EOF(3) THEN ' Get line color & store with x=0
  515.            INPUT #3, y%(i): LNi = LNi + 1: LINE INPUT #3, LN$(LNi)' Save line name
  516.            IF y = -1 THEN GOSUB LoadLabels ' All labels listed at end of file
  517.            END IF
  518.         LOOP: nmp = i  'nmp points to 0,-1 that ends all data (but the value
  519.                        'of X% and y% are 0,0 until file is saved.
  520.      LET CDY = LATcen: CDX = LONcen'Center display on ORIGIN
  521.      LET CPX = CDX: CPY = CDY 'Cursor Posn to Center of Display
  522.      LET Z = 2: LNptr = 1: REM start at first point and first line segment
  523.      CLOSE #3: RETURN:     REM first X% value is map color.  2nd val is 1st pt
  524.  
  525.  
  526.  
  527. LoadLabels: k = 0
  528.      DO WHILE NOT EOF(3)
  529.         k = k + 1: INPUT #3, ML$(k), MLa(k), MLo(k), MLr(k)
  530.      LOOP
  531.      IF MLa(k) = 0 OR MLo(k) = 0 THEN nml = k - 1 ELSE nml = k
  532.      RETURN
  533.        
  534. SaveMap: GOSUB BoxLine23
  535.    PRINT "Enter file name to save if other than "; MapFile$;
  536.    INPUT a$: IF a$ <> "" THEN MapFile$ = a$
  537.    F$ = MapFile$
  538.    GOSUB BoxLine23: PRINT "Saving map to file named "; F$; " ..."
  539.    OPEN F$ FOR OUTPUT AS #4
  540.    IF Fault = 70 THEN CLOSE #4: GOTO SaveMap
  541.    PRINT #4, LATo; ","; LATtext$
  542.    PRINT #4, LONo; ","; LONtext$
  543.    PRINT #4, ppdV; ","; VS$
  544.    PRINT #4, LATcen; ","; LATcen$
  545.    PRINT #4, LONcen; ","; LONcen$
  546.    PRINT #4, MapRng; ","; MapRng$
  547.    PRINT #4, MinRng; ","; MR$
  548.    PRINT #4, TextLine$
  549.    j = 1
  550.    FOR i = 1 TO nmp
  551.        IF x%(i) <> 0 THEN WRITE #4, x%(i), INT((y%(i) / Yfactr) + .5)
  552.        IF x%(i) = 0 AND i = nmp THEN PRINT #4, " 0,-1"
  553.        IF x%(i) = 0 AND i <> nmp THEN
  554.           PRINT #4, "0,0"
  555.           PRINT #4, y%(i); ","; LN$(j): j = j + 1
  556.        END IF
  557.    NEXT i
  558.    PRINT #4, "0,"; LN$(LNi)
  559.    FOR k = 1 TO nml
  560.        PRINT #4, ML$(k); ","; : WRITE #4, MLa(k), MLo(k), MLr(k)
  561.    NEXT k: CLOSE #4: LOCATE 24, 1:
  562.    Changed = 0
  563.    IF nmp > MaxNumPoints OR nml > MaxNumLABELS THEN
  564.       CLS : LOCATE 9, 29: PRINT "CAUTION!": PRINT : PRINT
  565.       IF nmp > MaxNumPoints THEN
  566.          PRINT "            The number of points,"; nmp; "is greater than"; MaxNumPoints
  567.       END IF
  568.       IF nml > MaxNumLABELS THEN
  569.          PRINT "            The number of LABELS,"; nml; "is greater than"; MaxNumLABELS
  570.       END IF
  571.       LOCATE 18, 12
  572.       PRINT " Therefore this map will not work with APRS (yet) "
  573.       LOCATE 23, 1: INPUT "HIT Enter to continue..."; a$
  574.    END IF: GOTO DrwMPaCur
  575.  
  576. CurDrwMap: CDX = CPX: CDY = CPY: GOTO DrawMap: REM Re-center at CURSOR location
  577.  
  578. DrwMPaCur: CPX = CDX: CPY = CDY: GOSUB DrawMap
  579.            REM After drawing map, Put cursor at center
  580.            RETURN
  581.  
  582. DrawMap: IF USGS THEN RETURN
  583.     Display$ = "MAP": COLOR 15, 0
  584.    'Draw to range scale RS and center display CDX and CDY
  585.    'Original Map was 40 pix-per-deg Horiz and 20 vert for 200 display
  586.    'Now ppdH and ppdV are variables.  The scaling factor KP is 1 for
  587.    'the original map scale.
  588.    DO WHILE RS < 320 / ppdV: RS = RS * 2: LOOP
  589.    IF RS > 8192 THEN RS = 8192
  590.    KP = 100 * 100 / (RS * ppdV)'This is to scale it down from the 120 maps
  591.    Sfac = 50 * 200 / RS        'Till 307 had been 100*120
  592.  
  593.    Lfac = COS(CDY / 57.296)
  594.    Hfac = (640 / 350) * (3 / 4) * Lfac
  595.    dx = ppdV * (LONo - CDX)
  596.    dy = ppdV * (LATo - CDY)
  597.         
  598.    CLS : LOCATE 1, 2: PRINT "Redrawing Map"
  599.    REM first put ORIGIN and map CENTER on the map
  600.    LINE (320 - KP * dx, Ycen - KP * dy)-(960 - KP * dx, Ycen - KP * dy), 14
  601.    LINE (320 - KP * dx, Ycen - KP * dy)-(320 - KP * dx, 3 * Ycen - KP * dy), 14
  602.    CMX = 320 + Sfac * (CDX - LONcen) * Hfac'new
  603.    CMY = Ycen + Sfac * (CDY - LATcen) * Yfactr
  604.    LINE (CMX - 27, CMY)-(CMX + 27, CMY), 14
  605.    LINE (CMX, CMY - 20)-(CMX, CMY + 20), 14
  606.    CIRCLE (CMX, CMY), 10, 14
  607.    CIRCLE (320 - KP * dx, Ycen - KP * dy), 12, 14
  608.    s = 0: GOSUB MapPoint: REM Redraw MapPoint
  609.    StrtPt = -1
  610. DP: FOR i = s TO nmp - 1
  611.        x = 320 + KP * (x%(i) - dx) * Hfac
  612.        y = Ycen + KP * (y%(i) - dy) * Yfactr
  613.        X1 = 320 + KP * (x%(i + 1) - dx) * Hfac
  614.        Y1 = Ycen + KP * (y%(i + 1) - dy) * Yfactr
  615.        IF StrtPt = -1 THEN CIRCLE (x, y), 3, 9: StrtPt = 0
  616.        IF x%(i + 1) <> 0 THEN
  617.           IF RdsOn OR LineColor <> 12 THEN LINE (x, y)-(X1, Y1), LineColor
  618.           IF i = Z THEN SavClr = LineColor
  619.        ELSE
  620.           REM LINE (x - 3, y - 3)-(x + 3, y + 3), 10, B: StrtPt = -1
  621.           CIRCLE (x, y), 5, 10: StrtPt = -1
  622.           LineColor = y%(i + 1): i = i + 1
  623.           IF Display$ = "SHOW" AND LineColor > 8 THEN LineColor = LineColor - 8
  624.        END IF
  625.    NEXT i
  626.    GOSUB Cursor
  627.    GOSUB ReDraw
  628.    REM MapPoint went here
  629.    GOSUB DrawHist: REM draw GPS history track
  630.    IF Display$ = "SHOW" THEN
  631.       GOSUB ShowMaps
  632.    ELSE
  633.       LOCATE 25, 1: PRINT "Use +/- to move MAPpoint.  N/P for Next/Previous Feature.  H for HELP!.";
  634.       LOCATE 1, 61
  635.       PRINT "POINTS"; nmp; "= "; INT((nmp / MaxNumPoints) * 100); "%";
  636.       LOCATE 2, 61
  637.       PRINT "LABELS "; nml; "= "; INT((nml / MaxNumLABELS) * 100); "%";
  638.       LOCATE 3, 61: PRINT "CENTER  "; MID$(STR$(LATcen), 2, 5);
  639.       LOCATE 3, 75: PRINT MID$(STR$(LONcen), 2, 5)
  640.       LOCATE 4, 61: PRINT "SCALE   (ppd)"; ppdV
  641.       LOCATE 5, 69: PRINT "Range"; LEFT$(STR$(MapRng), 5)
  642.    END IF
  643.  
  644. labels:
  645.    IF Labls THEN
  646.       FOR i = 1 TO nml ' Now plot labels on map
  647.       IF RS <= MLr(i) OR Key$ = "S" THEN
  648.          LET x = 320 + Sfac * (CDX - MLo(i)) * Hfac'new
  649.          LET y = Ycen + Sfac * (CDY - MLa(i)) * Yfactr
  650.          IF Tags AND y > 14 * Yfactr AND y < 350 * Yfactr AND x > 8 * (LEN(ML$(i)) + 1) AND x < 632 THEN
  651.             LOCATE y / (14 * Yfactr), (x / 8) - LEN(ML$(i)): PRINT ML$(i);
  652.             END IF
  653.          END IF
  654.       NEXT i
  655.    END IF
  656.    GOSUB ShowMap: RETURN
  657.          
  658. ShowMap: REM this shows the map boarder of the loaded map
  659.     x = 320 + KP * (CDX - LONcen) * ppdV * Hfac'new
  660.     y = Ycen + KP * (CDY - LATcen) * ppdV * Yfactr
  661.       by = MapRng * Sfac * Yfactr / 60
  662.       bx = by * 640 / (400 * Yfactr) * Lfac'old
  663.       C = 15
  664.     LINE (x - bx, y - by)-(x + bx, y + by), C, B
  665.     RETURN
  666.  
  667. Cursor: CIRCLE (CUX, CUY), 4, 0
  668.      CUX = 320 + Sfac * (CDX - CPX) * Hfac'new
  669.      CUY = Ycen + Sfac * (CDY - CPY) * Yfactr
  670.      CIRCLE (CUX, CUY), 4, 14
  671.      x = INT(CPX): y = INT(CPY): Xm = (CPX - x) * 60: Ym = (CPY - y) * 60
  672.      x$ = RIGHT$(STR$(x), 3) + " "
  673.      LOCATE 1, 2: PRINT "RNG"; RIGHT$("   " + STR$(RS), 4) + "      "
  674.      LOCATE 2, 2: PRINT "LAT "; y; MID$(STR$(Ym) + "   ", 2, 5)
  675.      LOCATE 3, 2: PRINT "LON "; x$; MID$(STR$(Xm) + "   ", 2, 5)
  676.     
  677.      LOCATE 24, 1: PRINT "Cursor coordnts:"; TAB(17);
  678.      PRINT INT(.5 + dx + (CUX - 320) / KP); TAB(23); INT(.5 + dy + (CUY - Ycen) / KP);
  679.      REM LOCATE 24, 55: PRINT "Degrees: ";
  680.      REM PRINT LEFT$(STR$(CPY) + " ", 7); LEFT$(STR$(CPX) + "   ", 7);
  681.      LOCATE 1, 16: PRINT "Decimal";
  682.      LOCATE 2, 15: PRINT LEFT$(STR$(CPY) + " ", 8);
  683.      LOCATE 3, 15: PRINT LEFT$(STR$(CPX) + "   ", 8);
  684.      LINE (0, 0)-(178, 42 * Yfactr), 12, B'Box around it
  685.      LINE (0, 0)-(116, 42 * Yfactr), 12, B'Box around it
  686.      LET a$ = "": LET B$ = "": RETURN
  687.  
  688. HELP: CLS : COLOR 15, 1: LINE (0, 0)-(639, 18 * Yfactr), 14, BF
  689.       LOCATE 1, 20: PRINT " MAPFIX.bas HELP SCREEN Ver "; Ver$
  690.       LOCATE 3, 1
  691.       
  692.       PRINT " The cursor is shown in LAT/LON, map offset and decimal degrees.   The ORIGIN,"
  693.       PRINT " CENTER and BORDER are shown (but only the CENTER and RANGE in MAPLIST.map are"
  694.       PRINT " actually used by APRS.  Labels are right justified to the point just after the"
  695.       PRINT " last letter.  CALLS & OBJECT names will be left justified."
  696.       PRINT ""
  697.       PRINT " OPERATIONS          MAP FUNCTIONS         @N - NEW FEATURE   LABEL COMMANDS"
  698.       PRINT " H - HELP SCREENS    @C- Change CENTER     @A - ADD point      S - SHOW labels"
  699.       PRINT " D - map DIRECTORY   ^C- Change COLOR      @D - DELETE point  @L - add a LABEL"
  700.       PRINT "^S - SAVE MAP!!!     M - MAPLIST.apr file  @K - Kill feature   L - LABELS off"
  701.       PRINT " R - RESET pointers  O - OTHER map bordrs  @M - MOVE point"
  702.       PRINT " Q - QUIT            @R- set map RANGE     @T - TRIM borders"
  703.       PRINT "                                               "
  704.       PRINT " MAP COMMANDS        POINTER MOVEMENTS   USGS CD ROM CMDS   DIGITIZER & GPS "
  705.       PRINT " SPACE to draw map   N - Next Feature    B - BOX PPD area   @O- OPEN COMMS"
  706.       PRINT " ARROWS  (shft)      P - Prev Feature    U - USGS overlay   @B- BEGIN new line"
  707.       PRINT " PgUP/DN (ctrl)      G - Go to Pointer   T - Test smoother "
  708.       PRINT " HOME to Cursor      F - Find point      @S- SMOOTH file"
  709.       PRINT " END to map center   +/- move POINTER    @U- USGS BUILD!    @G- GPS OVERLAY"
  710.       PRINT "                                         @J- JOIN lines"
  711.       PRINT "                                         ^R- REDRAW on/off"
  712.       PRINT : LINE (0, 190 * Yfactr)-(639, 190 * Yfactr), 15
  713.               
  714.       IF Display$ <> "HELP" THEN
  715.        LOCATE 25, 1
  716.        PRINT " HIT H AGAIN FOR MORE HELP SCREENS, OR SPACE BAR FOR MAP...";
  717.       END IF
  718.       Display$ = "HELP"
  719.       LINE (0, 0)-(634, 348 * Yfactr), 15, B
  720.       RETURN
  721.  
  722.  
  723.  
  724. REM ************* HERE IS THE CODE BROUGHT IN FROM APRS  ***************
  725.  
  726. LdMapLst: GOSUB BoxLine23: INPUT "FileSpec for MAPLIST.apr if not \APRS\MAPLIST.APR"; a$
  727.           IF a$ <> "" THEN F$ = a$ ELSE F$ = "\aprs\Maplist.apr"
  728.           OPEN F$ FOR INPUT AS #3: IF Fault <> 0 THEN RETURN
  729.           i = 1: NumGood = 0
  730.           INPUT #3, DfltY: LINE INPUT #3, a$
  731.           INPUT #3, DfltX: LINE INPUT #3, a$
  732.           INPUT #3, BestRng: LINE INPUT #3, a$: DfltR = BestRng
  733.           INPUT #3, GMToffset: LINE INPUT #3, a$
  734.     WHILE a$ <> "* BEGIN *": LINE INPUT #3, a$: WEND ' Skip comment block
  735.           REM RS = BestRng: REM center display
  736.           REM RS = 2 ^ INT(LOG(RS) / LOG(2))'Rng is intgr power of 2
  737.           REM CPX = CDX: CPY = CDY 'Cursor Posn to Center of Display
  738.     WHILE NOT EOF(3) AND i <= UBOUND(MapName$)
  739.           INPUT #3, MapName$(i), LATcen(i), LONcen(i), MapMax(i)
  740.           LINE INPUT #3, Comment$(i)' IGNORE ALL comment fields
  741.           REM now ignore maps that start with a *
  742.           IF LEFT$(MapName$(i), 1) <> "*" THEN NumGood = NumGood + 1
  743.           NumMaps = i: i = i + 1
  744.     WEND: CLOSE #3
  745.     IF NumGood >= MaxNumMAPS - 1 THEN
  746.           CLS : LOCATE 2, 5
  747.           PRINT "WARNING: Too many ACTIVE MAPS (more than"; MaxNumMAPS; ") in MAPLIST.map file for APRS"
  748.           LOCATE 4, 10: PRINT "Use EDITOR to suppress mapnames with an (*) that you don't need."
  749.           PRINT : PRINT : PRINT : MapListLoaded = -1
  750.           INPUT "HIT RETURN to continue"; a$
  751.     END IF
  752.     RETURN
  753.  
  754. ListMAPlist: IF NOT MapListLoaded THEN GOSUB LdMapLst
  755.    GOSUB ListHeader
  756.    FOR i = 1 TO NumMaps
  757.        IF i / 19 = INT(i / 19) THEN
  758.           LOCATE 25, 1: PRINT "HIT RETURN to continue"; : INPUT a$
  759.           GOSUB ListHeader
  760.        END IF
  761.        PRINT MapName$(i); TAB(14);
  762.        PRINT INT(LATcen(i) * 100) / 100; TAB(21); INT(LONcen(i) * 100) / 100;
  763.        PRINT TAB(29); MapMax(i); TAB(36); LEFT$(LTRIM$(Comment$(i)), 43)
  764.    NEXT i
  765.    
  766.    LOCATE 25, 1: PRINT "LIST COMPLETE. CONTINUE WITH NEXT MAPFIX COMMAND...";
  767.    RETURN
  768.  
  769. ListHeader: CLS
  770.    PRINT "MAPS in MAPLIST.map (*MAPS are suppressed)     [For now, use EDITOR to modify]"
  771.    PRINT :
  772.    PRINT "MAP NAME      LATcen LONcen  RANGE COmments"
  773.    PRINT "------------  ------ ------- ----- -------------------------------------------"
  774.    RETURN
  775.  
  776. DrwAndShow: IF NOT MapListLoaded THEN GOSUB LdMapLst
  777.             Display$ = "SHOW": GOSUB DrwMPaCur
  778.  
  779. ShowMaps: IF MapSize > RS / 2 THEN MapSize = RS / 2
  780.     LOCATE 25, 1: PRINT " Drawing all maps >"; MapSize;
  781.     PRINT "mi.  F3 to see smaller, F4 for bigger, SPACE to cancel.";
  782.     LINE (0, 336 * Yfactr)-(639, 349 * Yfactr), 14, B
  783.     FOR i = 1 TO NumMaps
  784.    
  785.     x = 320 + Sfac * (CDX - LONcen(i)) * Hfac
  786.     y = Ycen + Sfac * (CDY - LATcen(i)) * Yfactr
  787.       dy = MapMax(i) * Sfac * Yfactr / 60
  788.       dx = dy * 640 / (400 * Yfactr) * Lfac
  789.       C = 15
  790.       IF MapMax(i) > 32 THEN C = 14
  791.       IF MapMax(i) > 64 THEN C = 12
  792.       IF MapMax(i) > 128 THEN C = 11
  793.       IF MapMax(i) > 256 THEN C = 13
  794.           
  795.     IF MapMax(i) > MapSize THEN
  796.        LINE (x - dx, y - dy)-(x + dx, y + dy), C, B
  797.        IF y + dy > 14 * Yfactr AND y + dy < 350 * Yfactr THEN
  798.           IF x + dx > 8 * (LEN(MapName$(i)) + 1) AND x + dx < 632 THEN
  799.              LOCATE (y + dy) / (14 * Yfactr), (x + dx) / 8 - LEN(MapName$(i))
  800.              IF MapMax(i) > RS / 4 THEN PRINT MapName$(i);
  801.           END IF
  802.        END IF
  803.     END IF
  804.     NEXT i: RETURN
  805.         
  806. REM ******* here is the code added by W7KKE for overlyaying track histoiries
  807. 'This module retrieves GPS history files so you can check the accuracy of
  808. 'the map
  809.  
  810. Hstdir: CLS : PRINT "HISTORY FILES DIRECTORY": PRINT
  811.          PRINT "To display HST files, please enter the path to your xxxxxxx.HST files."
  812.          PRINT "For example, the default '\APRS\*.HST' will show all maps in the APRS"
  813.          PRINT "directory.  Similarly '*.hst' will search your present QBasic directory."
  814.          PRINT "For any other path, enter the full file specification.": PRINT
  815.          
  816.          PRINT "Enter Filespec for searching the DIRECTORY (\aprs\*.hst)";
  817.          INPUT F$: IF F$ = "" THEN F$ = "\aprs\*.hst"
  818.          IF INSTR(F$, ".") = 0 THEN F$ = F$ + ".HST"
  819.          PRINT : PRINT : FILES F$
  820.          RETURN
  821.  
  822. LoadHst: GOSUB BoxLine23
  823.        INPUT "Which history file to load (ENTER for list, Q to quit)"; F$
  824.        IF UCASE$(F$) = "Q" THEN RETURN
  825.        IF F$ = "" THEN GOSUB Hstdir: GOTO LoadHst
  826.        a = INSTR(3, F$, "."): IF a = 0 THEN F$ = F$ + ".hst"
  827.        Fault = 0: F$ = UCASE$(F$): OPEN F$ FOR INPUT AS #3
  828.    
  829.        IF Fault = 53 OR Fault = 62 THEN Fault = 0: RETURN
  830.        GOSUB BoxLine23: PRINT "Loading track history from "; F$
  831.        
  832.     DO WHILE NOT EOF(3)
  833.        i = i + 1
  834.        INPUT #3, a$
  835.        HLAT(i) = VAL(MID$(a$, 26, 2)) + (VAL(MID$(a$, 28, 5)) / 60)
  836.        HLONG(i) = VAL(MID$(a$, 35, 3)) + (VAL(MID$(a$, 38, 5)) / 60)
  837.        maxhist = i
  838.      LOOP
  839.      CLOSE #3: Histloaded = -1
  840.      GOSUB BoxLine23: PRINT "File loading is complete.  GPS data is plotted."
  841.      REM fall through...
  842.  
  843. DrawHist:  'put history track on map
  844.      IF Histloaded THEN
  845.         size = 3: IF RS < 2 THEN size = size * 2 / RS
  846.         FOR i = 1 TO maxhist
  847.             HMX = 320 + KP * (CDX - HLONG(i)) * ppdV * Hfac'new
  848.             HMY = Ycen + KP * (CDY - HLAT(i)) * ppdV * Yfactr
  849.             CIRCLE (HMX, HMY), size, 13
  850.         NEXT i
  851.      END IF
  852.      RETURN
  853.  
  854. NewMap: CLS : PRINT "BEGINNING A NEW MAP FROM SCRATCH...": PRINT
  855.    PRINT "All points in an APRS map are measured as an"
  856.    PRINT "offset to the right and down from an origin."
  857.    PRINT
  858.    INPUT "Enter the LATITUDE  of the ORIGIN in decimal degrees"; LATo
  859.    INPUT "Enter the LONGITUDE of the ORIGIN in decimal degrees"; LONo
  860.    PRINT
  861.    PRINT "Choose the number of pixels per degree to set the map scale:"
  862.    PRINT
  863.    PRINT "Approximate size    Range from center  Pixels/Deg"
  864.    PRINT "----------------    -----------------  ----------"
  865.    PRINT "Big state or region        250             120"
  866.    PRINT "Typical state              100             300"
  867.    PRINT "Several Counties            50             600"
  868.    PRINT "Typical VHF range           25            1200"
  869.    PRINT "City streets (7.5 min maps) 12            2400"
  870.    PRINT
  871.    INPUT "Enter desired Pixels/Deg"; ppdB
  872.    IF ppdB = 0 THEN GOTO NewMap
  873.    REM In following lines, 500 is half of 999 (maximum nominal value for pts)
  874.    LATcen = LATo - (500 * Yfactr / ppdB)
  875.    LONcen = LONo - (500 / ppdB)
  876.    GOSUB StartMap: ppdV = ppdB
  877.    CLS : PRINT "YOU ARE NOW READY TO DRAW A NEW MAP...": PRINT : PRINT
  878.    PRINT "A white border has been drawn around the maximum size permitted for this map"
  879.    PRINT
  880.    PRINT "USING CURSOR WITHOUT DIGITIZER:  Move coursor to starting point for a NEW"
  881.    PRINT "feature and hit ALT-N.  Then enter new feature name (for reference purposes)"
  882.    PRINT "and continue moving cursor to the next point and hit ALT-A to add more points."
  883.    PRINT "Continue in this fashion, using ALT-N whenever you want to begin a NEW feature."
  884.    PRINT
  885.    PRINT "USING A DIGITIZER:  First, use ALT-O once to OPEN the digitizer COM port.  Then"
  886.    PRINT "use ALT-B to BEGIN each new map feature.  Enter the name and color of the new"
  887.    PRINT "feature.  Then use the digitizer mouse to add more points."
  888.    PRINT : PRINT : PRINT
  889.    PRINT "Add LABELS on the map at the current cursor location by using the ALT-L key. "
  890.    PRINT
  891.    PRINT "When you are finished, be sure to SAVE the map using the CTRL-S command..."
  892.    PRINT : PRINT : PRINT
  893.    PRINT "FOR HELP, REMEMBER THE  H  KEY!"
  894.    PRINT : PRINT : PRINT "Hit ENTER to proceed..."; : INPUT a$
  895.    RETURN
  896.  
  897. StartMap: REM This called by NEW and in middle of USGS build
  898.    LATcen$ = "LAT of CENTER": LONcen$ = "LON of CENTER"
  899.    MapRng = 60 * 500 * Yfactr / ppdB: REM 500 is half of full map size
  900.    MapRng$ = "Map range from center"
  901.    VS$ = "Pixels per degree"
  902.    MinRng = 1: MR$ = "No longer used"
  903.    TextLine$ = "NEW Map generated by MAPFIX.bas routine..."
  904.    IF Key$ = "NEW" THEN RS = 2 ^ INT(LOG(MapRng) / LOG(2))'Rng is intgr of VERTrng
  905.    CDX = LONcen: CDY = LATcen: CPX = CDX: CPY = CDY
  906.    nmp = 1: nml = 0
  907.    LNi = 1: LN$(1) = "Labels begin here"
  908.    RETURN
  909.  
  910.  
  911. DigiInit: CLS : PRINT : Digitizer = -1
  912.    PRINT "This routine will replace many CURSOR functions with the Digitizer's MOUSE."
  913.    PRINT "Assuming your digitizer can output an X,Y,C format."
  914.    PRINT
  915.    PRINT "Only Mercator projection charts will give absolutely accurate results.  Other"
  916.    PRINT "types, Lambert Conformal, Conical, etc will induce distortions."
  917.    PRINT
  918.    PRINT "It has not been tested with East Longitude or South Latitude."
  919.    PRINT : PRINT
  920.    PRINT "The digitizr should operate at 9600,N,8,1 in POINT mode with 200 LPI resolution."
  921.    PRINT "The FORMAT outputs X,Y,C values separated by commas (C is for button pressed."
  922.    PRINT
  923.    PRINT "Set up the digitizer according to your model's instructions.  For the model"
  924.    PRINT "23360, use the drawing board menu by pressing the mouse button 0 on the SETUP"
  925.    PRINT "label so that the LED is ON.  Then move the mouse to each other label and"
  926.    PRINT "use the 0 button to toggle the value ON or off as follows:"
  927.    PRINT
  928.    PRINT "POINT is ON                             PARITY 7/8 and 1 are ON "
  929.    PRINT "BAUDRATE 3 is ON                        FORMAT is ON ON off ON"
  930.    PRINT "DATA RATE doesn't matter                RESOLUTION off off ON"
  931.    PRINT : PRINT
  932.    INPUT "Is DIGITIZER connected to COM1 or COM2 (1)"; a$
  933.    IF a$ <> "2" THEN a$ = "COM1" ELSE a$ = "COM2"
  934.  
  935.    Port$ = a$ + ":9600,N,8,1,cs0,ds0,cd0"
  936.    OPEN Port$ FOR RANDOM AS #1
  937.    
  938.  
  939.    CLS : PRINT "FIRST LETS TEST THE DIGITIZER, AND GET THE MAP ON STRAIGHT.": PRINT
  940.    PRINT "Move your mouse (or pen) and hit the 0 button (or touch tablet) to see if the"
  941.    PRINT "digitizer is outputting in the desired format.  While doing this, it is a good"
  942.    PRINT "idea to verify that your map is on straight.  The Y values from the mouse"
  943.    PRINT "should give the same values for the same LATITUDE line on both the right and"
  944.    PRINT "left edges of the map.  If not, move your map to get it horizontal."
  945.    PRINT
  946.    PRINT "OUTPUT FORMAT:"
  947.    PRINT
  948.    PRINT "XXXXX,YYYYY,APn (Only the X and Y values are used (4 or 5 digits is ok)"
  949.    PRINT
  950.    LOCATE 25, 1: PRINT "Hit ENTER and press 0 button on mouse to continue...";
  951.    LOCATE 13, 1
  952.    DO UNTIL INKEY$ <> "": LINE INPUT #1, a$: PRINT a$: LOOP
  953.   
  954.    CLS : PRINT
  955.    PRINT "NEXT YOU MUST ESTABLISH THE SCALE OF YOUR DIGITIZER."
  956.    PRINT
  957.    PRINT "The scale is established by two points, the first near the"
  958.    PRINT "upper left corner, the second near the lower right corner."
  959.    PRINT
  960.    PRINT "To get the best accuracy on maps not exactly MERCATOR, use points within the "
  961.    PRINT "area where you are working, not on the extreme corners.  IE:  choose points"
  962.    PRINT "that are in the center of the upper left quadrant and the lower right quadrant."
  963.    PRINT
  964.    PRINT "To establish the upper left reference point:"
  965.    INPUT "     Enter lat  (deg,min)"; LATref1, M: LATref1 = LATref1 + M / 60
  966.    INPUT "     Enter long (deg,min)"; LONref1, M: LONref1 = LONref1 + M / 60
  967.    PRINT
  968.  
  969.    PRINT "Place the mouse on the upper left point and press the 0 button."
  970.    LINE INPUT #1, a$: SOUND 150, 3
  971.         digix1 = 5000 - VAL(LEFT$(a$, 5))
  972.         digiy1 = VAL(MID$(a$, 7, 5))
  973.         PRINT "Digitizer reads "; digix1, digiy1; " for this point.": PRINT
  974.  
  975.    PRINT "NOW Establish the lower right reference point:"
  976.    INPUT "      Enter lat  (deg,min)"; LATref2, M: LATref2 = LATref2 + M / 60
  977.    INPUT "      Enter long (deg,min)"; LONref2, M: LONref2 = LONref2 + M / 60
  978.    PRINT
  979.    PRINT "Place digitizer pen on lower right point."
  980.    LINE INPUT #1, a$: SOUND 150, 3
  981.         digix2 = 5000 - VAL(LEFT$(a$, 5))
  982.         digiy2 = VAL(MID$(a$, 7, 5))
  983.         PRINT "Digitizer reads "; digix2, digiy2; " for this point.": PRINT
  984.  
  985.    REM Find delta lat/long between reference points
  986.    REM Calculate degrees per x/y unit
  987.         degx# = (LONref1 - LONref2) / (digix1 - digix2)
  988.         degy# = (LATref1 - LATref2) / (digiy1 - digiy2)
  989.  
  990.    CLS : PRINT "YOU ARE NOW READY TO USE THE DIGITIZER TO ENTER POINTS INTO MAPFIX..."
  991.    PRINT
  992.    PRINT "The digitizer works just about like the cursor and arrow keys in MAPFIX.  Any"
  993.    PRINT "point identified by the digitizer will be ADDED just as if you had hit ALT-A."
  994.    PRINT "All points are added to a feature after the current MapPoint identified by the"
  995.    PRINT "white circle. "
  996.    PRINT
  997.    PRINT "With the digitizer, do NOT use the ALT-N NEW command which always begins at the"
  998.    PRINT "current cursor location.  For the digitizer, use ALT-B to BEGIN a new feature."
  999.    PRINT "You will be asked to identify the name and color of the new feature.  From then"
  1000.    PRINT "on, just move the digitizer mouse (or pen) to ADD new points.  "
  1001.    PRINT
  1002.    PRINT "If your digitizer mouse has 4 buttons, use the first (left) button for ADDing "
  1003.    PRINT "points, use the 4th (right) button to just move the cursor with no action."
  1004.    PRINT : PRINT
  1005.    INPUT "Hit ENTER to continue with MAPFIX..."; a$
  1006.    RETURN
  1007.  
  1008. GetXY: LINE INPUT #1, a$: SOUND 150, 3
  1009.        a = INSTR(a$, ","): IF a = 0 THEN RETURN
  1010.            x = 5000 - VAL(LEFT$(a$, a - 1))
  1011.        B = INSTR(a + 1, a$, ","): IF B = 0 THEN B = LEN(a$)
  1012.        y = VAL(MID$(a$, a + 1, B - (a)))
  1013.        Btn = VAL(RIGHT$(a$, 1))
  1014.        CPY = ((y - digiy2) * degy#) + LATref2
  1015.        CPX = ((x - digix2) * degx#) + LONref2
  1016.        IF LOC(1) <> 0 THEN a$ = INPUT$(LOC(1), #1)'Clear input buffer
  1017.        RETURN
  1018.  
  1019. BoxPPD: GOSUB BoxLine23: INPUT "Enter the desired PPD"; a$: ppdB = VAL(a$)
  1020.    dy = (30000 / ppdB) * Sfac * Yfactr / 60
  1021.    dx = dy * 640 / (400 * Yfactr) * Lfac
  1022.    GOSUB BoxLine23: LINE (CUX - dx, CUY - dy)-(CUX + dx, CUY + dy), 13, B
  1023.    PRINT "The box represents the largest APRS map that can be made with that scale."
  1024.    RETURN
  1025.         
  1026.  
  1027. GetUSGS: REM This used for both U=OVERLAY and by ALT-U = USGS BUILD!
  1028.    ReDraw = 0: USGS = -1: ni = 0: nt = 0: j = 0: NumLines = 0: LE = 1: OE = 1
  1029.    IF Key$ <> "U" THEN
  1030.       IF ppdB <> 0 THEN ppdV = ppdB
  1031.       LATo = CDY + (500 * Yfactr / ppdV)
  1032.       LONo = CDX + (500 / ppdV)
  1033.       GOSUB BoxLine23: PRINT "Improve LAT ORIGIN of "; LATo; : INPUT LATo
  1034.       GOSUB BoxLine23: PRINT "Improve LON ORIGIN of "; LONo; : INPUT LONo
  1035.       GOSUB BoxLine23: INPUT "LATitude extent (100%)"; a$
  1036.             IF a$ <> "" THEN LE = VAL(a$) / 100
  1037.       GOSUB BoxLine23: INPUT "LONgitude extent (100%)"; a$
  1038.             IF a$ <> "" THEN OE = VAL(a$) / 100
  1039.       dx = ppdV * (LONo - CDX)
  1040.       dy = ppdV * (LATo - CDY)
  1041.       KP = 100 * 100 / (RS * ppdV)
  1042.       LATcen = CDY: LONcen = CDX: GOSUB StartMap
  1043.       LATtext$ = "Decimal LATITUDE  of map ORIGIN"
  1044.       LONtext$ = "Decimal LONGITUDE of map ORIGIN"
  1045.       TextLine$ = "Map generated by MAPFIX from USGS 2,000,000:1 CD ROM (data valid mid-1980's)"
  1046.    END IF
  1047.    Lmax = 500 + 600 * LE: Lmin = 500 - 600 * LE 'Max=1100 and Min =-100
  1048.    Omax = 500 + 600 * OE: Omin = 500 - 600 * OE
  1049.    GOSUB BoxLine23: INPUT "Enter path and filename of XTRACTED USGS data file"; a$
  1050.    a = INSTR(a$, "."): IF a = 0 THEN a$ = a$ + ".GRF"
  1051.    OPEN a$ FOR INPUT AS #3
  1052.    IF Fault <> 0 THEN RETURN
  1053.    GOSUB BoxLine23: INPUT "Enter RANK cutoff. (all roads = 99) for ST and WB use 20, or 10"; a$
  1054.    MaxRnk = VAL(a$)
  1055.    REM PRINT "raw data format.....", "   LineID", "#-Rnk-Atbts", "  NumPts"
  1056.    LOCATE 5, 67: PRINT "RANGE: "; INT(30000 / ppdV): LOCATE 24, 1
  1057.    IF Key$ = "U" THEN
  1058.       PRINT "While USGS OVERLAYED, do not redraw map or you will have to do it again...";
  1059.    ELSE PRINT "Blue circles start lines, Green Box ends.  Red points discarded, Yellow Kept!";
  1060.    END IF
  1061.  DO UNTIL EOF(3) OR LNi = MaxNumLines - 1
  1062.    NumLines = NumLines + 1
  1063.    LOCATE 1, 61: PRINT "TOTAL POINTS: "; nt
  1064.    LOCATE 2, 61: PRINT "POINTS USED:  "; ni
  1065.    LOCATE 3, 61: PRINT "TOTAL LINEs:  "; NumLines
  1066.    LOCATE 4, 61: PRINT "LINEs USED:   "; LNi
  1067.    a$ = INPUT$(20, 3): REM PRINT a$;
  1068.         LnID$ = LEFT$(a$, 7)
  1069.         Rank$ = MID$(a$, 8, 2): Rank = VAL(Rank$): LOCATE 6, 67: PRINT "RANK: "; Rank
  1070.         Npts$ = MID$(a$, 10, 6): Npts = VAL(Npts$)
  1071.         AtCd$ = MID$(a$, 16, 5)
  1072.         a$ = LTRIM$(LnID$) + "-" + Rank$ + "-" + AtCd$
  1073.         REM PRINT , LnID$, a$, Npts$
  1074.    IF Rank < 24 THEN SavClr = 4 ELSE SavClr = 7
  1075.    IF Rank < 20 THEN SavClr = 12
  1076.    IF Rank < 14 THEN SavClr = 10
  1077.    LineOK = 0: IF Key$ <> "U" THEN GOSUB BeginF
  1078.    FOR i = 1 TO Npts
  1079.     a$ = INPUT$(20, 3): IF Rank > MaxRnk THEN GOTO Skp
  1080.     REM IF VAL(Rank$) > 99 THEN GOTO Skp
  1081.     LA = VAL(LEFT$(a$, 2)) + VAL(MID$(a$, 3, 2)) / 60 + VAL(MID$(a$, 5, 2)) / 3600
  1082.     LO = VAL(MID$(a$, 8, 3)) + VAL(MID$(a$, 11, 2)) / 60 + VAL(MID$(a$, 13, 2)) / 3600
  1083.     IF Key$ = "U" THEN
  1084.         REM Following lines used to limit points if just doing an OVERLAY only
  1085.         IF LA > CDY + RS / 60 OR LA < CDY - RS / 50 THEN GOTO Skp 'off screen
  1086.         IF LO > CDX + RS / 40 OR LO < CDX - RS / 40 THEN GOTO Skp
  1087.         REM s$ = MID$(a$, 16, 5)
  1088.         REM PRINT S$, LA, LO
  1089.         END IF
  1090.     y% = (LATo - LA) * ppdV: x% = (LONo - LO) * ppdV
  1091.     IF Key$ <> "U" AND (x% > Omax OR x% < Omin) THEN GOTO Skp'this ignores points off PPD
  1092.     IF Key$ <> "U" AND (y% > Lmax OR y% < Lmin) THEN GOTO Skp'scale
  1093.     LineOK = -1
  1094.        X1 = 320 + KP * (x% - dx) * Hfac
  1095.        Y1 = Ycen + KP * (y% - dy) * Yfactr
  1096.        IF i > 2 THEN
  1097.           REM LINE (x, y)-(X1, Y1), 6
  1098.           dd = LO - LOb: IF dd = 0 THEN dd = .0000001
  1099.           dn = LA - LAb
  1100.           s = dn / dd' Note that 1>s>.01 for Xdelta of 1 to 100
  1101.           IF ABS(s) < .1 AND ABS(Lsp) < .1 THEN
  1102.              sd = 1
  1103.           ELSEIF ABS(s) > 10 AND ABS(Lsp) > 10 THEN sd = 1
  1104.           ELSEIF ABS(dd) < .004 AND ABS(dn) < .004 THEN sd = 1
  1105.           ELSEIF s <> 0 THEN sd = Lsp / s
  1106.           ELSE sd = 0
  1107.           END IF
  1108.           IF sd > 2 OR sd < .5 OR i = Npts THEN
  1109.              CIRCLE (x, y), 1, 14
  1110.              IF Key$ <> "U" THEN GOSUB KeepLine
  1111.           ELSE CIRCLE (x, y), 1, 4
  1112.           END IF
  1113.           Lsp = s: nt = nt + 1
  1114.        ELSE Lsp = 0: CIRCLE (X1, Y1), 4, 9
  1115.           IF Key$ <> "U" THEN GOSUB KeepLine 'keeps first two lines
  1116.        END IF
  1117.        LAb = LA: LOb = LO
  1118.        x = X1: y = Y1
  1119.        
  1120. Skp: NEXT i
  1121.      IF Key$ <> "U" THEN
  1122.         IF LineOK THEN nmp = nmp - 1: Z = Z - 1: ni = ni + 1 ELSE GOSUB CanclF
  1123.      END IF
  1124.      LINE (x - 3, y - 3)-(x + 3, y + 3), 10, B ' Last Point
  1125.  LOOP
  1126.  IF LNi > MaxNumLines - 2 THEN LOCATE 12, 20: PRINT "PROCESSING STOPPED... TOO MANY LLINES!..."
  1127.    CLOSE #3
  1128.    RETURN
  1129.  
  1130. KeepLine: x%(Z) = x%: y%(Z) = y%: nmp = nmp + 1: Z = Z + 1: ni = ni + 1: RETURN
  1131.  
  1132.  
  1133. Scrunch: i = 0: Pt = 0: nt = 0: ni = 0: GOSUB BoxLine23
  1134.     INPUT "Enter slope filter ratio 1.2 to 5 (typically 1.5)"; a$
  1135.     IF a$ = "" THEN slope = 1.5 ELSE slope = VAL(a$)
  1136.     DO UNTIL i >= nmp - 1
  1137.        i = i + 1
  1138.        X1 = 320 + KP * (x%(i) - dx) * Hfac
  1139.        Y1 = Ycen + KP * (y%(i) - dy) * Yfactr
  1140.        IF x%(i) <> 0 THEN
  1141.           Pt = Pt + 1
  1142.           IF Pt > 2 THEN
  1143.              LINE (x, y)-(X1, Y1), 6
  1144.              dd = x - X1: IF dd = 0 THEN dd = .01
  1145.              dn = y - Y1
  1146.              dst = ((dd * dd) + (dn * dn)) ^ .5
  1147.              s = dn / dd' Note that 1>s>.01 for Xdelta of 1 to 100
  1148.              IF s = 0 THEN s = .05
  1149.              IF ABS(s) < .2 THEN s = .2 * SGN(s)
  1150.              IF ABS(s) > 5 THEN s = 5 * SGN(s)
  1151.              IF ABS(s) <= .2 AND ABS(Lsp) <= .2 THEN
  1152.                 sd = 1
  1153.              ELSEIF ABS(s) >= 5 AND ABS(Lsp) >= 5 THEN sd = 1
  1154.              ELSE sd = Lsp / s
  1155.              END IF
  1156.              IF ABS(dd) > 50 * KP OR ABS(dn) > 30 * KP THEN sd = 0
  1157.              REM IF ABS(dd) < 5 OR ABS(dn) < 4 THEN sd = 1
  1158.              IF sd > slope OR sd < 1 / slope OR x%(i + 1) = 0 OR NumRej > 4 THEN
  1159.                 ni = ni + 1: CIRCLE (x, y), 2, 14: NumRej = 0
  1160.                 
  1161.              ELSE CIRCLE (x, y), 1, 4: NumRej = NumRej + 1
  1162.                 IF Key$ <> "T" THEN
  1163.                    i = i - 1: nmp = nmp - 1
  1164.                    FOR ii = i TO nmp
  1165.                        x%(ii) = x%(ii + 1): y%(ii) = y%(ii + 1)
  1166.                    NEXT ii
  1167.                 END IF
  1168.              END IF
  1169.              Lsp = s: nt = nt + 1
  1170.           ELSE Lsp = 0: nt = nt + 1: ni = ni + 1: CIRCLE (X1, Y1), 4, 9
  1171.           END IF
  1172.        ELSE Pt = 0: nt = nt + 1: ni = ni + 1
  1173.           LOCATE 1, 61: PRINT "TOTAL POINTS: "; nt
  1174.           LOCATE 2, 61: PRINT "SAVED POINTS: "; ni
  1175.        END IF
  1176.        x = X1: y = Y1
  1177.     LOOP
  1178.     RETURN
  1179.  
  1180. Join: REM Search for end=begin point values and CONCATONATE if equal!
  1181.    LNptr = 0: i = 0: k = 0: GOSUB BoxLine23: PRINT "Lines joined: ";
  1182.    DO UNTIL i >= nmp
  1183.       i = i + 1
  1184.        IF x%(i) = x%(i + 2) AND y%(i) = y%(i + 2) AND y%(i + 1) = LColor THEN
  1185.           nmp = nmp - 2: LNi = LNi - 1: k = k + 1: LOCATE 23, 15: PRINT k
  1186.           FOR j = i + 1 TO nmp: x%(j) = x%(j + 2): y%(j) = y%(j + 2): NEXT j
  1187.           FOR j = LNptr TO LNi: LN$(j) = LN$(j + 1): NEXT j
  1188.        ELSEIF x%(i) = 0 THEN LColor = y%(i): LNptr = LNptr + 1
  1189.        END IF
  1190.    LOOP: GOTO DrawMap
  1191.  
  1192. END
  1193.  
  1194.